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1 Introduction 

The phenomenon of self-organized criticality is widely studied in Physics 
from different perspectives. (We refer toP,[2], [H], [22], [9], [15], [IT], [12], 
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[5] [I3|, [ID] , [IS] , [E] for various studies). Roughly speaking it is the property 
of systems to have a critical point as attractor. 

The purpose of this paper is to analyze this phenomenon in the frame 
work of stochastic evolution equations. To the best of our knowledge this is 
the first time this is done in the presence of a stochastic force and in such 
generality in a mathematically strict way. Let us introduce our framework. 

Let G be an open bounded domain of M"^, (i = 1, 2, 3, with smooth bound- 
ary dG . We shall study here the nonlinear stochastic diffusion equation, 

' dX{t) - A<il{X{t))dt 3 a{X{t))dW{t), in (0, oo) x G, 
< ^(X(t))9 0, on (0, oo) X (9^, (1.1) 

X{0,x) = X on G, 

where x is an initial datum and : M — 2^^ is a maximal monotone (possibly 
multivalued) graph with polynomial growth and random forcing term 

oo 

a{X)dW = f^kXdpk Ck, t > 0, 
fc=i 

where {cfc} is an orthonormal basis in L^(^), {/ifc} is a sequence of positive 
numbers and {(3^} a. sequence of independent standard Brownian motions on 
a filtered probability space (fi, {^t\t>Oi P)- 
We note that (t{X) is defined by 

oo 

a{X)h = Y,^^kX{Keu)2ek, ^ h e L'{G), 

k=l 

where (■, ■)2 is the scalar product in L^(^). 

The equation models the dynamics of flows in porous media and more 
generally the phase transition (including melting and solidification processes) 
in the presence of a random forcing term a{X)dW. 

Existence for stochastic equations of the form (11.11) with additive and 
multiplicative noise was studied in [6] under the main assumption that is 
monotonically increasing, continuous and such that 

f ^(0) = 0, *'(r) < ai|rp~i + 02, VreM, 

(1.2) 

^(s)ds > a3|rP+^ + a4, V r G M, 
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where ai > 0, 03 > 0, 02, 04 > and m > 1. (See also [7] and [2T] for general 
growth conditions on 

Here we shall study equation (11.11) under the following assumptions. 

Hypothesis 1.1 (i) is a maximal monotone multivalued function from 
R into R such that G ^(0). 

(a) There exist C > and m > 1 such that 

snp{\9\ : 6 e *(r)} < C(l + |r|™), V r G M. 

(Hi) The sequence {fik} is such that 

00 

fc=i 

where Afc are t/ie eigenvalues of the Laplace operator —A in G with 
Dirichlet boundary conditions. 

We recall that the domain of Ais H^{^)n Hl{ff). 
Since for x G H~^{ff) 

\xek\\ < Ci|efc|^2(^) \x\\ < Ci\l\x\\ (1.3) 

and hence 

00 00 

\W{x)\\l2{L^{0),H-H0)) = ^AWkW <Ci^lll\l\xW, (1.4) 

k=l k=l 

it follows by (lii) that a{x) G L2{L'^{ff), H-^{ff)) (the space of all Hilbert- 
Schmidt operators from L?{iff) into H~^{ff)) and that it is Lipschitz con- 
tinuous from H~^{0') into L2{L^{^), H^^{iP')). Under these assumptions we 
shall prove that if x G U'{ff),p > max{2m, 4}, then there is a unique strong 
solution to equation (11. ip which is nonnegative if so is the initial data x. With 
respect to the situation considered in [7] , in the present case one does not as- 
sume that the range of \l/ is all of R, which is quite unusual for porous media 
equations. Also Hypothesis ll.ll (^i) allows monotonically increasing functions 
\1/ with a finite number of discontinuities (jumps), ri,...,r7v. One must, of 
course, fill the jumps by taking \l/(rj) = [^{vj + 0), \l/(rj — 0)], j = Vi, ...,riy. 
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It should be mentioned that several physical problems with free boundary 
and with phase transition can be put into this functional setting. For instance 
if 

{ai{x — a), for X < a 
[0,p], forx = a (1.5) 
a2{x — a) + p, for x > a, 

with a, p, tti, 0^2 £ (0, +oo), then f 1 1.1 1) models the phase transition in porous 
media or in heat conduction (Stefan problem). If \&(x) = p sign x where 
p > and 

f ^, ifx^O 
sign X = < (1.6) 

[ [-1,1], ifx = 0, 

then (11. ip reduces to the nonlinear singular diffusion equation 

dX{t) - p div {6{X{t))VX{t))dt = a{X{t))dW{t), 

where 6 is the Dirac measure concentrated at the origin. 
Other examples such as the Heavside step function 

f 0, if X < 
H{x) = < [0,1], if X = 
( 1, if X > 0, 

or \l/(x) = 1x1°" sign x with < a < 1 also satisfy Hypothesis II. 1[ 
In particular the equation 

dX{t) - A{H + X){X{t) - Xc)dt = a{X{t) - Xc)dW{t), (1.7) 

where A > 0, represents the continuous, stochastic version of the Bak, Jang, 
Wiesenfeld sand pile model [1]. (See [1] for a deterministic presentation of the 
model.) This is a diffusion problem with free boundary driven by a random 
forcing term proportional to X{t) — Xc where Xc is the critical density and 
X{t) is the density at the moment t. 

Taking into account the numerical simulation in 1-D (see [2]), one might 
expect that the time evolution of the system displays self-organized critical- 
ity, i.e. the supercritical region {X{t) > Xc} is absorbed asymptotically in 
time by the critical one {X{t) = Xc}. Here we shall prove that, e.g. in all 
examples fll.5l) - fll.7p above, this indeed takes places with high probability 
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under appropriate assumptions on the parameters and more precisely that 
the supercritical region "vanishes" into the critical one in finite time with 
high probability, at least if fi^ = for all > + 1 for some iV e N. We 
emphasize that this is in particular true when the noise is zero. In this case 
one gets an explicit bound for the time when this happens (cf. Remark 14.41 
below) . 

The plan of this paper is the following. The main results are presented 
in Section 2 and are proven in Section 3. In Section 4 we prove a finite time 
extinction type result for solutions to (11. ip which displays a self-organized 
criticality behaviour. 

The following notations will be used. L^^ff), p > 1, is the usual space 
of p-integrable functions with norm denoted by | ■ \p. The scalar product 
in L^((^) and the duality induced by the pivot space L^(^) will be denoted 
by (■, ■)2- H^{ff) C L^(^), k = 1,2, are the standard Sobolev spaces on 
iff, while Hl{ff) is the subspace of H^{0') with zero trace on the boundary. 
For p,g G [l,+oo] by LlyliO^T); U>{n; H)) {H a Hilbert space) we shall 
denote the space of all g-integrable processes u : [0, T] — > LP{fl] H) which are 
adapted to the filtration {^t}t>o- 

By Cvi/([0, T]; L^(f2; H)) we shall denote the space of all if- valued adapted 
processes which are mean square continuous. L{H) denotes the space of 
bounded linear operators equipped with the usual norm. 

In the following by H we shall denote the distribution space 

H = H~\ff) = {Hl{ff))' 
endowed with the scalar product and norm defined by 

where A = -A with D{A) = H\ff) n Hl{ff). 

In terms of A equation (II. ip can be formally rewritten as 

dX{t) + A^{X{t))dt 3 a{X{t))dW{t), 

(1.8) 

X{0,x) = X. 

Its exact meaning will be precised later (see Definition 12. II below) . 
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It should be recalled, however, that the operator x — A^{x) with the 
domain 

{x e L^{ff) n H'^{ff) : there is 7] e ifo(^),r/ G ^(x) a.e. in ff} 

is maximal monotone in H := H~^{ff) (see e.g. [3]) and so the distribution 
space H offers the natural functional setting for the porous media equa- 
tion (11. ip or its abstract form (11.81) . However, the general existence theory 
of infinite dimensional stochastic equations in Hilbert space with nonlinear 
maximal monotone operators (see [TT], [20]) is not applicable in the present 
case and so a direct approach must be used. 

Fnally, in this paper we use the same letter C for several different positive 
constants arising in chains of estimates. 

2 Existence, uniqueness and positivity 

Definition 2.1 Let x ^ H . An H-valued continuous ^t-adapted process 
X = X(t,x) is called a solution to (11.11) {equivalently (11.81) ) on [0,T] if 

X G LP{n X (0, T) X ^) n L\0, T; L\n, H)), p> m, 

and there exists rj G ^^/'"(fi x (0,T) x ff) such that P-a.s. 

(X(t,x),ej)2 = {x,ej)-2+ / / r]{s,^)Aej{^)d^ds 

Jo Jff 

oo „t 

+ y^fik {X{s,x)ek,ej)2d/3k{s), Vj gN, tG [0,r], 
k=i 

(2.1) 

rj G ^{X) a.e. in x (0,T) x ^. (2.2) 

Below for simplicity we often write X{t) instead of X{t^x). 

From the stochastic point of view the solution X given by Definition 12.11 
is a strong one, but from the PDE point of view it is a solution in the sense 
of distributions since the boundary condition \E'(X)^0 on d&' is satisfied in 
a weak sense only. 

Theorem 12.21 below is the main existence result. 
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Theorem 2.2 Assume that d = 1,2,3 and that Hypothesis holds. Then 
for each x G U'{0'), p > max{2m, 4} there is a unique solution X G 
L^(0, T; ^)) to (11.11) . Moreover, if x is nonnegative a.e. in G then 

P-a.s. 

X{t,x){i) > 0, for a.e. e (0, oo) x G. 

As mentioned earlier. Theorem 12.21 was proven in [0] for a differentiable 
\1/ satisfying conditions (11.21) and for p > max{m + 1,4}. It should be said, 
however, that in contrast with what happens for coercive functions \i/ arising 
in [6], here it seems no longer possible to extend the existence result to all 
X G H-^{ff), x>0. 



3 Proof of Theorem 2.2 



We shall consider the approximating equation 

dXx{t) + A{^x{Xx{t)) + XXx{t))dt = a{Xx{t))dW{t), 

(3.1) 

Xx{0,x) = X, 
where A > and 

^x{x) = i (x - (1 + A^)-i(x)) G ^((1 + A^)-i(x)) 
A 

is the Yosida approximation of We recall that \E'a is Lipschitzian and 
monotonically increasing and so x — > "^xix) + Ax is strictly monotonically 
increasing and bounded by Ci(l + |x|™) and {"^xix) + Ax)x > A|xp for all 
X G M. By P, Theorem 2.2] (applied with m = 1), for each x G H~^(G) 
equation (13. ip has a unique solution 

X;, G L\Q X (0,T) X G)nLl,{Q,C{[0,T];H)) 

in the sense of Definition 12.11 Here as usual C{[0,T]; H) is equipped with 
the supremum norm. Moreover, ( see e.g. [20l Theorem 4.2.5]) the following 
Ito formula holds 



E\Xx{t)\l, + 2E / f {^x{Xx{s)) + XXx{s))Xx{s)d^ds 
Jo J ff 



X 



oo „t 

.i + J^/i^E / \Xx{s)ek\\ds. 
k=i 



(3.2) 



We note that since 

\Xxek\-i < C|efc|//2(^')|XA|_i < CXk\Xx\-i, 
(cf. ([IS])) we have by Hypothesis O^iii) (cf. 

^/i^E / \Xx{s)euWds<C¥. \Xx{s)Wds. (3.3) 
fc=i 

Lemma 3.1 There exists a constant C > such that for all p > 2 and all 

X G 



ess.sup,,[o,T]E|X;,(t,x)|^<exp(^C^^ |a;|^, V A > 0. (3.4) 



(3.5) 



(3.6) 



Proof. We know from |6l Lemma 3.4] (with m = 1) that as e ^ 
XI ^ Xx strongly m L^(0, T; L\n- H)), 

XI Xx weakly star m L^(0, T; LP(fi; LP{ff))), 
where X| is the solution to the approximating equation 

' dXiit) + {Ax),Xi{t)dt = a{Xi{t))dW{t), t > 0, 
Xm = X, 

where 

Axx = A{-^x{x) + Xx) = -A{-^x{x) + Ax), 
D{Ax) = {xeHn L\^) : ^x{x) + Xx e Hi 
and {Ax)e is the Yosida approximation of Ax, 

{Ax)e = \{I-{I + eAxr^), £>0. 

Furthermore, by [6l Lemma 3.2] we have that X| e L'^{Vt]C{[Q,T];L'^{ 
As a matter of fact the results of [6] were proven for smooth nonlinear func- 
tions while \E'a is only Lipschitz; but the extension to lipschitzian functions 
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\1/ satisfying f 1 1.2 1) is immediate. In fact, one might take a smoother ap- 
proximation of for instance the moUifier * p\ (pa(^) = \ p{^/'^)iP ^ 
C^(M),p > 0, / pdr = 1) which still remains monotonically increasing and 
has all properties of \E'a. 

Next we apply Ito's formula (13. 6p for the function (p{x) = ^ \x\^. More 
precisely, we first apply Ito's formula to v^-y(x) = ^ |(1 + 7 > 0, 

and then we let 7^0. We have (for details see the proof in [6], Lemma 3.5]), 

E^(XKt)) +E !\{A,),Xl{s), \Xl{s)r'Xl{s)),ds 
Jo 

= cpix) + / \Xlis)r'\Xlis)e,\'d^ ds d^ (3.7) 

k=i 

< ^{x) + ^ CE f [ \Xl{s)\Pd^ ds, 
^ Jo Je 

since by Sobolev embedding |efc|oo < C\k for all G N. If is the solution 
to the equation 

- eH^x{Y^) + \Yl) = XI vi/,(y,^) + XY^ e Hl{ff), 

then (see [HI (3.25)]) \Y^\p < \Xf\p and therefore 

{{A,),Xl \Xir'Xl), = \{Xi- Yl \Xir'Xl), > 0. 
Then by (13.71) it follows, via Gronwall's lemma, that 

nm)\i<\Ai exp (c^), 

where C is independent of x, A and t. Now one obtains (13.41) by letting e 
tend to and taking into account (13.51) . □ 

From now on let us assume that p > max{4, 2m} and x G U'{G'). From 
Lemma [3. II it follows that for a subsequence {A} — we have 

' Xx-^X weakly in LP{Vt x (0, T) x ^), 

and weakly star in in L~(0, T; Lp(1]; Lp(^))), 

' *a(Xa) ^ 7] weakly in LP/™(fi x (0,T) x ^), ^"^'^^ 

in particular in /.^(fi x (0,T) x ^), 
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because by Hypothesis (zi), 

|^A(a;)| < < C(l + |xr), VxeM. 

is the minimal section of By fl3.4p we have for A — ^ 

XXx strongly in L^in x (0, T) x (3.9) 
Clearly X and rj are adapted processes. On the other hand, we have 
d{Xxit) - X^{t)) - A($,(X,(t)) - ^,iX^{t)) + XX^it) - fiX^{t))dt 

= (a(X,(t)) - a{X^mdW{t) 

and therefore once again applying Ito's formula (cf. (13.21) ) we obtain for 
a>0,te [0,T], 

l\X,{t)-X,m-ie-''' 



+ 








+ (AXa(s) - iiX^{s)){Xx[s) - X^(s)) 



e-^^'di ds 



<{C^filXl-i^a] I \Xxis-X,{s))We-''' ds + MxAt), VA,/i>0, 

(3.10) 
where 

MxAt) ■= f e~^'{Xx{s) - X^{s),a{Xx{s) - X^{s))dW{s))2 
Jo 

is a real local valued martingale. To derive (I3.10p we used that x = X'^xlx) + 
(1 + A^)"^(x) and thus for all x,y eR 

{^^I,(x)-<^I,{y)){x-y) = [M>A(x)-vI/^(y)][(l + AM/)"i(x)-(l + ^M/)"i(y)] 

+ [vl>,(x)-v[/^(y)][Avl>,(x)-/iVl>^(y)], 
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and that the first summand on the right hand side is nonnegative because \l/ 
is monotonically increasing and "^xix) e ^'((l + A\l/)~^(x)). Hence for a > 
large enough we obtain for all A, /i G (0, 1) and t E [0, T] 

l\X,{t)-X,m\e--' 

<Cmax{A,/i} f I (|vl>,(X,(s))p + + |vl>^(X^(s))|2 (3.11) 



Hence by the Burkholder-Davis-Gundy inequality (for = 1) we get for all 
A,/ie(0,l),re[0,T], 

sup |X,(t)-X^(t))|2^e-"* 



<Cmax{A,^}E / / (|^a(Xa(s))P + |Xa(s)|2 + |^^(X^(s))p (3.12) 
+ |X^(s)|2)e-"Xrfs + CE^VA(s)-X^(s)|iie-2"^ds^ ^ . 



But 

1/2 

E( / |XA(s)-X^(s)|V-^°^ci/ 



\ 1/2 



<E sup |X;,(s)-X^(s))|_ie-t-f r|XA(s)-X^(s)|V""'^^s') 

se[0,r] VJo / 

< ^ E sup |Xa(s) - X^(s))|^e-"^ + CE T |X,(s) - X^(s)|^e-"^cis. 

4 sG[0,r] Jo 

(3.13) 

Taking into account that by Hypothesis ll.l( ii) 

|*a(Xa)| <C(l + |XAr), VA>0, 
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and that by ([S3D {^a} is bounded in (fi x (0, T) x G) for p > max{4, 2m}, 
we infer by f l3.12p . (13.131) and Gronwall's lemma that {^a} is a Cauchy net 
in L2(n;C([0,T];iJ)) Hence for A ^ 

Xx-^X inL2(fi;C([0,T];i/)). (3.14) 

In order to complete the proof of the existence part of Theorem 12.21 it suffices 
to show that 

r/(cj,t,0 e a.e in X (0,T) X ^. (3.15) 

Since the operator 

LP(fix(0,T)x^) ^ L™(fix(0,T)x^) c (l]x (0, T) x ^), X ^(X), 
in the duahty pair 

'vi^i X (o,r) X ff),LP{Q X (o,r) X ff)' = L^in x (o,t) x e 

is maximal monotone, it suffices to show that (see e.g. [3]) 

hminfE I I ^!x{X>)Xxdidt<¥. [ [ riXd^dt. (3.16) 
Jo Jff Jo Jff 

To prove (13.161) we first note that by (13.21) we have 



liminfE / / ^xiXx)Xxd^dt+-E\X{t)\'i^ 
Jo 2 

1 1 °° /■* 

fc=i 



(3.17) 



because by (II. 3p . |(Xa — X)efc|_i < C\k\Xx — X\_i and so by Hypothesis 
0(iii) 

lim V/x^E / |XA(s)efc|^irfs = VyU^E / |X(s)efc|^ids. 
Next letting A tend to zero in (13.11) and using (13.81) we see tliat P-a.s., for all 

t e [0,T], 

ft °o f 

{X{t),ej)2 = {x,ej)2+ {ri{s),Aej)2ds + ^i2k {X{s)ek, ej)2d(3k{s). 

Jo f^^-^ Jo 

(3.18) 
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Note that by continuity the P-zero set does not depend on t G [0,T], since 
^/ifc / {X{s)ek,e,)^df3k{s)= / {e,,a{X{s))dW{s))^. 

k=l "^0 "^0 

In order to get (13.181) we have used the fact that by fl3.14p we have 
E 



{Xx{s)ek,ej)2d(3k{s)ds - / {X{s)ek,ej)2d(3k{s)ds 

Jo 

E [ {{Xx{s) - X{s))ek, e,)lds < CXjXlT\Xx - X\ 
Jo 



2 

L2{Q,C{[0,T];H)) 



and therefore 



hmV'/ife / {Xx{s)ek,ej)2d/3kds = y2l^k {X{s)ek, ej)2d/3kds. 
^^\=i ^0 Jo 

Therefore (13.181) follows and this yields, via Ito's formula (applied to (X(t), Cj)!, 
t G [0,T]) and summation over j that 



-E\X{t)W+E f [ TjXdids 
2 Jo J e 

1 1 /•* 

= -E|x|^ + - ^/i^E |X(.)e,|^rf., VtG [0,T]. 



(3.19) 



Comparing (I3.17P and (13.191) we get (I3.16P . Hence X is a solution to (II. ip as 
claimed. 

To prove uniquenss we take two solutions X'^^'^ and X*^^) with correspond- 
ing rj^^^ and 77*^^^. Repeating the argument above we obtain 

^ E\X^^\t) - X^^\t)W 



+E / / (r/(i)(s) - r]^^\s)){X^^\s) - X^^\s))did 



s 



JO 



-^Y.^l^ \ \{^^^\s)-^^^\s)yk\\ds, VtG[0,T]. 
fc=i 
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Since, because \l/ is monotone, the second term on the left is positive, by 
(11.31) . Hpothesis ll.l( iii) this imphes X^^^ = X^'^^ by Gronwall's lemma. 

Finally, if x > a.e. in we know by [6l Theorem 2.2] that Xx > 
P-a.s. and so by flSTTll) it follows that X >0, a.e infix (0, T) x ^ as desired. 
This completes the proof of Theorem 12.21 □ 



Remark 3.2 Theorem 12.21 extends to any dimension d > 1 if one modifies 
condition (Hi) in Hypothesis 1 1 . 1 1 as in [HI Condition 4.1]. 

Remark 3.3 The existence part of Theorem 12.21 remains true for stochastic 
porous media equations with additive noise, i.e. 



dX - A<i/{X)dt = VQ dW{t), 
where \1/ satisfies Hypothesis 11.11 and 

oo 
k=l 

with 

oo 
fc=l 

The proof is exactly the same and so, it will be omitted. 

Proposition 3.4 Let Xx,X E (0,1), be as above, x E L^{^). Assume that 
satisfies Hypothesis with m = 1 and for some 6 > 0, 

{x-y){x-y)>6ix-yf, V (x, x), (y, G ^. (3.20) 

ThenXx,X E L^^i^.T- L'^iVt, Hl{ff))) and 

limE|X,-X|i,(o_^^^.(^))=0. (3.21) 

Proof. A simple calculation reveals that 

{^xix) - ^x{y)){x -y)>-\x- y\\ V x, y e M 
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for A sufficiently small. Then defined hj'^xir) := x{r) — ^ r, r G M, is 
increasing and so by Ito's formula we have 

E\X,it)\l + ^ \X,{s)\lr^^^/s < a (3.22) 

As a matter of fact, we shall apply Ito's formula not directly to equation (13.11) 
but to equation (13. 6p (cf. the proof of Lemma [3.11 to obtain ( 13. 7p ). Thus we 
get 

^ E\Xml + Ej\{A,%Xiis),Xiis))2ds <l\x\l + CeJ^ \Xl{s)\lds. 
Next we have 

{{Ax)sXi, Xih = {Axil + eAx)-'Xl (1 + eAxr'Xl)^ + e\{Ax%Xi\l. 

Taking into account that Ax = A(\E'a + A/) and that r ^'a(^) — Sr/2 is 
monotonically increasing we get 

{iAx)eXlXlh > ^ ljVil + eAxr'Xl\'d^ + e\iAx)sXi\l 

Hence ^ 

E f \il+eAx)-'Xlis)\l,^^^^ds<C 

J 

and letting e — > we get (I3.22p and the first assertion (taking also into 
account (13. 5p ). 

To prove the second part we note that 

d{Xx - X^) - Al^xiXx) - ^^{X^) + XXx -fiX^ + ^6 {Xx - X^,)]dt 

= {aiXx)~aiX,))dW. 
Hence exactly the same arguments to derive (13. lip lead to 

i \Xx{t) - XM-ie-"-' + ^ £ \Xx{s) - X,{s)\le-^ds 
< Cmax{A,^} 1^ (\^x{Xx{s))\l + \^,{X,{s))\l 

+ \Xx{s)\l + \X,{s)\iy-"'ds + MxAt), 
15 



for a large enough and \, jj, & (0,1), t G [0,T]. Since m = 1, we have 
|^a(3^)| < ^"(1 + for all X G M, A G (0, 1), hence taking expectation we 
get 

^e/ |X,(s)-X^(s)|2rfs<Cmax{A,/i}E / + |X^(s)nrfs. 

By Lemma [3.11 with p = 2 and (13.81) this implies fl3.2ip . □ 

Besides Hypothesis 11.11 we shall now assume the following 

(iv) \E'(r) = p sign r + \i/(r), /or r G M, where p > 0, \E' : R — > M zs 
Lipschitzian, G C^(]R\{0}) and for some 6 > it satisfies \&'(r) > (5 
/or allr G M \ {0}. 

Here the signum is defined by (11.61) . 

Below we shall use an approximation to which is slightly different from 
defined before. Namely, below we consider 

^A(r) := p (sign)A(r) + ^(r), r G M, 

where (sign)A is the Yosida approximation of the sign, i.e. 

r 1 if r > A 
(sign)A(r) := } ^ if r G [-A, A] 
[ -1 if r < -A. 

We shall use the symbol \E'a also for this approximation and denote also by X\ 
the corresponding solution of (13. ip . This approximation in the special case 
of condition (iv) is much more convenient. We emphasize that all previous 
results remain true for this modified approximation. The proofs are the same 
and some parts even simplify. We therefore shall use all previous results for 
^^A and Xx as above without further notice. 

Proposition 3.5 The solutions Xx to (13. ip and X to (II. ip satisfy all con- 
ditions of Proposition \3.4\ and in addition 



E 



r [ \V (sign) x{Xx)Mdt < C, V A > 

Jo J £> 



and consequently rj G -^^^(0, T; L^(f2; 
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Proof. We set ^ 

9xir) := / [sign) x{s)ds, r G R, 
Jo 

and choose ip\ G C^(]R) such that 

(i) ^a(O) = 0. 

(ii) f'x{r) = J for |r| < A, (p'x{r) = 1 + A for r > 2A, (p'x{r) = — 1 — A for 
r < -2A. 

(hi) < (^^(r) < f for all r G M. 

It is easily seen that such a function exists and can be constructed simply by 
smoothing the function (sign) a- Let us denote the resulting function by fx. 
Then define 



pr 

/ fx{s)ds, r G 
Jo 



As mentioned above the arguments of the previous proofs extends to the 
present situation in order to prove that {Xx} is convergent to the solution 
X to (O). 

Now we shall apply Ito's formula to equation (13.11) (or, more exactly, to 
(13.61) and then let e ^ as in the proof of Proposition 13. 4p with "^x defined 
as above and to the function J^ipx{Xx)dC,. 

Arguing as in the proof of Lemma [3. II to obtain (13.71) . we get (recall that 
Xx{t) G 

E / ^x{Xx{t))d^-E [ {A{sign)x{Xx{s)) + A^{Xx{s)),ip'x{Xx{s)))2 ds 
J & Jo 

)di + CY,^A^ / / ifl{Xx{s))\Xx{s)ek\^dids 
e ^^.=1 "'0 Je 



< 



Pt P 

ipx{x)d^ + AXC V/i^^,E / / lx{s,0\ek\'d^ds, 
e Jo Je 



where 1a is the characteristic function of the set {(s,^) : < |Xa(s,^)| < 
2A}. 
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Concerning the first line we note that, since ip'^^ and ^ are monotonically 
increasing while as seen earlier ^a(^) ^ -^o (^)! we have by the Green formula 
that 

(A$(Xa),(pU^a))2 = - / ^\X,)^l{X^)\VXy^\^di<Q. 

J e 

This yields 

E / / (V(sign)A(XA), V^'x{Xx))2 d^ds < C, V A G (0, 1). 
Jo Jff 

Taking into account that 

-(A(sign)A(XA),<^',(XA))2 = (V(sign)A(XA), V<^',(Xa))2 > 0, a.e. 
and that Vip'^{Xx) = { VXx on : |^a(s,OI < A} we get 

E [ [ \V{sign)x{Xx)\'d^ds <C, V A G (0, 1), 
Jo J 

because V(sign)A(XA) = ]y{Xx) if \Xx)\ < A and V(sign)A(XA) = if 
|Xa)| > A. 

Then we get the desired estimate and since also by fl3.22p 

E /" [ \\/^{Xx)\^d^ds <C, VAg(0,1) 
Jo J ff 

and (sign)A(XA) + ^{Xx) r] weakly in /.^(fi x (0,r) x ^) as A ^ we 
infer that r] G L^(0,T; L^{n; H^{^)) as claimed. □ 



4 Extinction in finite time and self-organized 
criticality 

In this section we shall prove a finite extinction property for solutions of 
(11.11) in 1-D for a special density dependent diffusion coefficient function 
However, Lemma [4.11 below can be proved without restriction on dimension. 
So, for the moment we remain in our general framework. 
For simplicity we choose the Wiener process 

N 

W{t) = J2f^kekPk{t), t>0, (4.1) 

k=l 
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where iV G N. 

Besides Hypothesis 11.11 we shall assume Hypothesis (iv) (see page 16), 

i.e. 

(iv) \E'(r) = p sign r + \l'(r), for r G M, where p > 0, : M. ^ M. is 
Lipschitzian, \Ef G C^(]R\{0}) and for some 6 > it satisfies \i/'(r) > S 
for allr G M \ {0}. 

Here the signum is defined by (11.61) . 
Now let T be the stopping time 

r = inf{t > : |X(t,x)|_i = 0}, 

where X(t,x),t > 0, is the solution to (11.11) given by Theorem 12.21 for x G 
LP(^), p > max{4,2m}. 

Lemma 4.1 Under assumptions (i)-(iv) we have 

X{t,x) = 0, fort > r, P-a.s.. 

Proof. Set A = -A, D{A) = H^{ff) n Hl{ff). Define p : [0,T] x 1] ^ 
C|(^;M)by 

N 

M^) := -5^/^fcefc/?fc(t), tG[0,r], 
fc=i 

and p : [0,T] ^ Cl^ff-^) by 

AT 

E2 2 

k=l 

Define 

F(t) = e^(*)X(t), t > 0. 

Let D{A) be equipped with the graph norm of A and let D{Ay be its dual 
space, hence 

D{A) C ifo(^) C L^{ff) C i7"^(^) C D{A)'. (4.2) 

It is easy to see that for all G ^2, t G [0, T] the function e^*-*''^-* is a multiplier 
both in D{A) and in if, hence e'^^^'^^^Az G -D(^)' is well defined for all 
z G L2(^) and F(t) G ii. 
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Claim. We have 



Y{t) = x+ / e 



1 r 

'^"^'^ Ar]{s)ds - - / fiY{s)ds, te[0,T], (4.3) 







where the fist integral on the right hand side is a Bochner integral in D{A)', 
the second by (13. 8p is one in U'{0') C L^(^). In particular a posteriori the 
first integal is in if, continuous in if as a function of t G [0, T], P-a.s. 

Proof of the Claim. Let ip G D{A). As before we shall use (■, ■)2 also for 
the extended dualizations with pivot space L'^{^) as the ones in (14. 21) . Then 



for t G [0,T] 



oo 



Furthermore, we have by Ito's formula for all ^ G ^ 




Now fix j G N. Then by the stochastic Fubini Theorem 
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By Ito's product rule and fl3.18p we hence obtain 



-J' 



(e„e'^(^V)2 (Ae„r/(s))2 ds 







N 



fc=l 



2 Jo 



fc=i 

AT 



^f^k {ej,X{s))2 (ej,efee^('V)2 



fc=l 

After summing over j G M the two stochastic terms cancel and the claim 
follows since ip G D{A) was arbitrary. 

Below we work for P-a.s. u & Q, u fixed. Hence all constants C appearing 
below may depend on u. 

Consider the solution Xx G L^(0, T; L2(fi, to equation IKT\\ . 

By Proposition 13.41 we have 

liniE|XA - X\l2(^Q^T.LHf?)) = 



and ^a(^a) e L'^{0,T; L'^{n, H^{ff))) because is Lipschitz. 
On the other hand,we have as in fl4.3p for Yx = e^Xx 

= e'^W Ar/A(t) - \ mYxit), V t > (4.4) 

where 

Vx{t) = ^xiXxit)) G H', 

It follows by (I32ID that 



limE|FA-l^lL(o,T;L^(^))=0 (4.5) 
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and therefore for some sequence A„ — 

lim \Yx„ - >^|l2(o,T;L2(^)) = a.e. on Q. (4.6) 

n-^oo 

Below we simple write A instead of A.„. Next we have by fl4.4p that 

=(r/.(t),A(e'^Wn(t))>^-i(/i(tm(t),n(t))2 a.e.tG[0,T]. 



2 



dt 

(4.7) 

Also we have (for simplicity we take p = 1) 
(r/,(t),A(e'^Wn(t)))2 

= ((sign)A (e-'^Wn(t)) + $(e-'^Wn(t)), A{e^^'^Yx{t))), 
= - I (V(sign), (e-'^Wn(t)),V(e'^Wn(t)))rfe 



vl>'(e"^WFA(t))(V(e-^WFA(t)),V(e^WFA(t)))rfe 

(ivn(t)r-iiA(t)r ivMi)r)iA(t,o^e 



vi.'(e-^wn(t))(ivrA(t)r - in(t)r iv/i(t)r)^ie, 

because for y G -ffg' 



r 0, on{y ^ (-A,A)}, 
V (sign)A {y) = < 

[ i Vy, on{ye(-A,A)}. 

(Here 1a is the characteristic function of {(^,t) G ^x[0,T] : |e-''^*'^)yA(t, 0)1 < 
A} and (■, ■) is the euclidean scalar product in M".) Since > 5 and 
^' G L°°(M), /i G C([0,T] X ^) this yields 

(r/A(t), A(e^WrA(t)))2 < C (1^(^)1^ + A) . (4.8) 
Hence (14.71) and Gronwall's lemma imply 

\Yxm < e^(*-^) {\Y,{s)\l + CAT) a.e. t > s. 
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Now taking into account 04.61) and letting A ^ we get 



Y{t)\l<e^^'~^^\Y{s)\l 



a.e. t > s. 



(4.9) 



If y(-) is L2(^) -continuous then (US]) holds for all s,t G [0,T], t > s. Taking 
in dini) s = r A T we get Y{t) = for alH > r A T and since T > was 
arbitrary for all t > r as claimed. So, we have to prove that Y is L^(^)- 
continuous on [0,T]. For this we recall that by Proposition 13.51 we have 



Then by equation (gJD we have ^ G L^{0,T; H-\^)) and so, since Y e 
L^(0, T; iJQ(^)) P-a.s. by Proposition 13. 4[ by a well known interpolation 
result (see e.g. |3]), we conclude that Y G C([0, T]; L^((^)). This concludes 
the proof of Lemma 14. 1[ □ 

For proving our extinction result we need ^ C M, i.e. d = 1. To be more 
specific let <^ = (0,7r). Then efc(^) = sin/c^, ^ G [0,7r], = k'^ and 
L^(0,7r) C H continuously, so 



r] e L\0,T;H^{ff)), P-a.s.. 



(4.10) 




(4.11) 



Theorem 4.2 Let x G L^(0,7r), p > max{2m, 4}, be such that 



where 



N 




(4.12) 



k=l 



Then, for each ri G N 




(4.13) 



where by Lemma\4-1\ we have 



t{u) = sup{t > : |X(t,x)|„i > 0}. 
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Proof. By condition (iv) we see that 

r^'(r)>p|r|, V r G M. (4.14) 

Consider the solution Xx G L'^{0,T; L'^^Q; H^{0,7i))) to equation fIXT]) . 
Then by first applying Krylov-Rozovskii's Ito formula (cf.[T9l Theorem 1.3.1] 
or e.g. [20l Theorem 4.2.5]) and then the classical Ito formula to the real 
valued semi-martingale \Xx(t)\'^i,t G [0,T], and the function 

(^^(r) = (r + ^2)^/^ rGM, 

we find 

^ 1 ^ , |X,(t)e,|^(|X,(t)|^ + e') - |(X,(t)e.,X,(t)),,n 
2^^' (|X,(t)|^+.^)3/2 

+ (a(X,(t))ciW^(t),^^(|X,(t)p_,)X,(t)) 

2 ^ {\Xx{t)\±i + e'r' 
- ^^ (|X,(tl|%+"e^)V2 + 2(a(^A(t))cil^(t), ^;(|X,(t)|^)X,(t)). 

(4.15) 

Here Cat is given by fl4.12p and 

N 

a{Xxit))dW{t) = Y,f^kXx{t)ekdPk{t). 

k=l 

Integrating over t and letting A ^ we see that the right hand side of (14.151) 
converges to the right hand side of (14.161) below in L'^[Q] C{[0, T]; H)). But 
by (IS3D,(IS3D, flXT^ . (Km and by Proposition the same is true for the 
left hand side with limit 

M\Xit)\\)-M\^\-l)+ f I ^}'\ V{s)dids. 



JA\X{s)W + eY/^ 
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Taking into account 02.21) and fl4.14p we altogether obtain 



{\X{t)t, + e^y/ 

- ^^ (|x(!)l%+"e2)i/2 + 2(a(X(t))rfl^(t), ^',{\X{t)t,)X{t)) . 
Consequently by Lemma 14.11 for alH > 

/•Mr 

+2 / MX(s))diy(s),<^;(|X(s)|^)X(s)), P-a.s., 
Jo 

where 7 is defined by (14.41) . 
Clearly, we have 

j•t^T \X{s)\ 
lim / ' 12 ds = t^T, p-a.s.. 

e-oJo + £2)1/2 

Now, letting e tend to zero we get 



|X(t)|„i + 7p(t A r) < + / \X{s)Uds 



+ / l[o,.](s)(a(X(s))ciiy(s),X(s)|X(s)|:}) P-a.s. 
Hence by a standard comparison result 



+ / e^-(*-^)l[o,](.)(a(X(.))dW^(.),X(.)|X(.)|:i). 
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(4.17) 



Taking expectation and multiplying by (p7) we obtain 

[ e-^^^P(r > s)ds < 
Jo PI 

Writing P(r > s) = 1 — P(r < s) we deduce that 



(r < t) > 1 - &i ( / e-^^'^'ds 
PI \Jo 



and flCTD follows. □ 

In particular Theorem 14.21 applies to self-organized criticality stochastic 
models ()1.7p 



' dX{t) - A(p sign (X(t) - + - Xc))dt 



N 



3 cr{X{t) - Xc) ^ fikekdPk, t>0, 



k=l 



(4.18) 



p sign {X{t) - Xc) + *(X(t) - Xc) 3 0, on 9[0, tt], 

I, X(0,a:) = X. 

Here the function \1/ is as in assumption (iv) and Xc G M. 
Corollary 4.3 Assume that 

\X - Xc\-l < plC^^, 
where Cn is as in 04.121) and 7 as in 04.111) . Then for each G N 

I ^ I —1 



P(t-c <n)>l 



PI 



e-^'^'ds 



(4.19) 



where 



Tc = inf{t > : \X{t) - Xc\^i = 0} = sup{t > : \X{t) - Xc\-i > 0}. 
and X = X{t,x) is the solution to 04.18P in the sense of Definition 2.1. 



26 



We note that equation fll.7l) reduces to (14.181) by shifting the Heavside 
function with Xc- 

One must notice that if a; > Xc, i.e. if the initial state is in the supercritical 
region then by positivity result in Theorem 12.21 we have X{t) > Xc, P-a.s. 
for all t > 0. This means that the state remains in the supercritical-critical 
region for all the time. However, by (14.191) if ^'^^^'^ is small, it reaches the 
critical state Xc with high probability in a finite time i.e. the supercritical- 
critical region is completely absorbed by the critical one in a finite time. 



Remark 4.4 Let us consider the deterministic case. Then (14.131) implies 
that r < if n > M^- so r < But this, of course, also follows directly 

from (I4.17p . since we assume Cn = 0. 
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